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Dynamic renormalization group (RG) methods were originally used by Forster, Nelson and 
Stephen (FNS) to study the large-scale behaviour of randomly-stirred, incompressible fluids gov- 
erned by the Navier-Stokes equations. Similar calculations using a variety of methods have been 
performed since, but have led to a discrepancy in results. In this paper, we carefully re-examine 
in d-dimensions the approaches used to calculate the renormalized viscosity increment and, by in- 
cluding an additional constraint which is neglected in many procedures, conclude that the original 
result of FNS is correct. By explicitly using step functions to control the domain of integration, 
we calculate a non-zero correction caused by boundary terms which cannot be ignored. We then 
go on to analyze how the noise renormalization, absent in many approaches, contributes an 0(k 2 ) 
correction to the force autocorrelation and show conditions for this to be taken as a renormalization 
of the noise coefficient. Following this, we discuss the applicability of this RG procedure to the 
calculation of the inertial range properties of fluid turbulence. 
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1. INTRODUCTION 

The large scale behaviour of randomly stirred fluids was originally studied by Forster, Nelson and Stephen (FNS) 
[l|, They used a dynamic renormalization procedure to explore the effects of the progressive removal of small 
(length) scales in a perturbative model under several types of forcing. As they note, their study is only valid at the 
smallest momentum scales, and as such the study is well below the momentum scale of the inertial range [3[. Later, the 
procedure used by FNS was extended by Yakhot & Orszag (YO) [H [1[ to a more general forcing spectrum (of which the 
studies of FNS were special cases) and used to calculate the energy spectrum and a value for the Kolmogorov constant 
in the inertial region. While their arguments allowing them to calculate inertial range properties are contested > 
these issues are not the main focus of this paper. Instead, we will concentrate on another disagreement related to the 
results for the renormalized viscosity and noise. 

In the papers of FNS and YO @, the authors calculate the viscosity increment, quantifying the effect of the 
removed subgrid scales on the super-grid scales. They find the prefactor Ad(e) = (d 2 — d — e) / 2d(d + 2) (from .5], with 
FNS in agreement for their specific cases of study) . The disagreement is centred around the use of a certain change 
of variables employed by FNS and YO. This substitution has been highlighted as a cause for concern (for example, 
since naively the symmetric domain of integration appears to be shifted, violating conditions for the identities 
used to be valid. Using methods that do not introduce any substitution, again for a general forcing spectrum, Wang 
k Wu (WW) [9] and Teodorovich [6j arrive at a different, incompatible result for the viscosity increment. Instead, 
they find the prefactor A* d = (d — l)/2(d + 2) = Ad(0). This e-free result is also used in the more field-theoretic 
work of Adzhemyan et al. [10]. Later, Nandy attempted to determine which of the results was correct using a 
"symmetrization argument" and agreed with the original (general forcing) result by YO. 

The method used by FNS and YO has found wide-ranging application, for example in soft matte r sy stems, such 
as the KPZ and Burgers' equations [T2rll5j , and the coupled equations of magnetohydrodynamics [l6|. Given the 
extensive use of this approach, it is unsatisfactory to have any lingering disagreement on the basic methodology. The 
aim of this paper is, therefore, to settle this dispute once and for all. There cannot be two different results for the 
same quantity. We will show that an extra constraint mentioned by FNS causes the elimination band not to be 
shifted, and that for substitution- free methods there are neglected boundary terms. These are evaluated and shown 
to compensate exactly the difference between A * d and A4. We then show how correct treatment does not require a 
symmetrization to obtain the Yakhot-Orszag result. 
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FIG. 1: The energy spectrum for a turbulent flow. Renormalization from A <C kd introduced by FNS in their IR study of 
randomly stirred flows takes you to the fixed point at k = 0. Iterative averaging from k c ~ O.lkd takes you to the non-Gaussian 
fixed point k* , which marks the beginning of a line of fixed points along fc _5//3 . 

In addition to renormalization of the viscosity, there is also renormalization of the noise. All treatments consider an 
input noise that is Gaussian with the forcing spectrum parametrized as Wok~ v , where k is the wavenumber associated 
with the force. At one-loop order each of the two vertices will have a factor of the inflowing momentum, thus leading 
to a k 2 contribution to the forcing spectrum. Both FNS and YO acknowledge this k 2 correction. In FNS they treat 
two specific cases, y = —2 and y = 0. In the former, they find a renormalization to Wq whereas, in the latter, 
they conclude that all higher order corrections are subleading. YO restrict their analysis to y > —2 and once again 
conclude all higher order corrections are subleading. We explicitly show how the leading contribution will always go 
as k 2 and as such can only be taken as an multiplicative renormalization for the case y = —2, as noted in [17j . We 
find the prefactor agrees with Ad(e) found by FNS and YO (with y = —2) and show it to be incompatible with the 
e-independent A* d . 

Another author, Ronis [l8j , calculates the viscosity and noise renormalization using a field-theoretic approach. His 
analysis agrees with FNS and YO for y = —2, although appears to be presented for general y. As we will argue in 
Sec. |H this seems unjustified as the noise is only renormalized for the case y = — 2. 

The paper is organised as follows. In Sec. [5J we give a brief discussion on the validity and limitation of this type of 
low-fc renormalization scheme for a fluid system. In Sec. [3] our calculation for viscosity renormalization is done and 
then in Sec. [4] for noise renormalization, along with comparison with other analyses. The results are summarised in 
table HI Finally, in Sec. [5J we present our conclusions along with a brief discussion of the relevance of this type of 
renormalization scheme for calculating inertial range quantities. 

2. DISCUSSION AND RELEVANCE OF APPROACH 

We start with a brief discussion on the region of validity of this method and its limitations. Turbulence is often 
viewed as an energy cascade, where energy enters large length scales in the production range and is progressively 
transferred to smaller and smaller scales, until viscous effects dominate and it is dissipated as heat. There must be a 
balance between the energy dissipated and the energy transferred through the intermediate scales, otherwise energy 
would build up and the turbulence would not remain statistically steady. Thus the dissipation rate, s, controls how 
small the smallest length scales need to be to successfully remove the energy passed down, giving the Kolmogorov 

scale r\ — [v z je) 1 ^ = 1/kd- When the Reynolds number is sufficiently large, there exists a range of intermediate 
scales where the energy flux entering a particular length scale from ones larger than it is the same as that leaving it 
to smaller ones and is thus not dependent on the wavenumber. This is the inertial range. 

The energy spectrum and a summary of the various ranges of it are presented in Fig. [T] (based on a similar figure 
in Q). In the RG approach, the smallest length scales (largest wavenumber scales k) are removed and an effective 
theory is obtained from the remaining scales. There is high- and low-energy asymptotic freedom since the renormalized 
coupling becomes weak in both limits Q. The dynamic RG method used by FNS introduces a momentum cutoff 
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A well below the dissipation momentum scale, below the inertial range even (see figure [I]), in the production range 
and removes momentum scales towards k — 0. As such, this method can only ever account for the behaviour on the 
largest length scales. The production range is highly dependent on the method of energy input, and so it is obvious 
that the properties of the lowest modes will also share this dependence. Taking the forcing to be Gaussian then allows 
Gaussian perturbation theory to be used, since the lowest order is simply the response to this forcing. Since the 
inertial range is highly non-Gaussian, we do not expect to study the inertial region with this analysis. 

An alternative RG scheme called iterative averaging (McComb [1, instead takes a cutoff k c ~ 0.1 kd and 

removes successive shells of wavenumbers down to a non-Gaussian fixed point k* , which marks the beginning of a 
line of fixed points following fc~ 5 / 3 through the inertial region (see figure [lj. The asymptotic nature of this method 
therefore cannot tell us anything about the forcing spectrum, and is only dependent on the rate at which energy is 
given to the system. No assumptions about Gaussian behaviour are made. 

Using the energy spectrum (figure [TJ, we see the location of the IR procedure of FNS/YO and how it is inapplicable 
for the calculation of inertial range statistics. Put simply, it does not have access to the inertial range, just as iterative 
averaging does not have access to the production range. This is discussed further in the conclusions, Sec. O 



3. CALCULATION 



The motion of an incompressible Newtonian fluid in d-spatial dimensions, subject to stochastic forcing, /, is 
governed by the Navier-Stokes equation (NSE) which, in configuration space, is 

du a (x,t) ( ( _ ^ d \ t„ +\_ ldp(x,t) | 2 



;jl y U ^ X,t ^dx^) Ua ^ X ' t ^ = ~'^~dx' ^ u oVu a (x,t) + f a (x,t), (1) 

where u(x,t) is the velocity field, p(x,t) is the pressure field, p is the density of the fluid and vq is the kinematic 
viscosity. The index a S {1, . . . , d} and there is an implied summation over repeated indices. We consider an isotropic, 
homogeneous fluid and, using the Fourier transform defined by 

the NSE may be expressed in Fourier-space as 

f d d j d£l f d d p dQ' 

(iu + v k 2 )u a (k 7 uj) = f a {k,u) + \ a M aM {k) J ^ d+1 J ^y+i fyuy(p, Q') 

(2ir) d+1 6{j +p- k)8(n + n' - uj) 



k a u a (k,uj) = , (4) 

where the incompressibility condition (V • u = 0) has been used to solve for the pressure field in terms of the velocity 
field. In equation (|3]), we have also introduced Ao (= 1) as a book-keeping parameter to the non-linear term and the 
vertex and projection operators, respectively, are defined as 



M a/3l (k) = — k/3P ai {k) + fc 7 P Q(3 (fc) 

k k 

P (k) - - a 1 

1 oryV, "v u ct7 jo ? 



(5) 



and contain the contribution from the pressure field. The integral over p, fi' could be trivially done to follow FNS 
0, YO [| and Wang & Wu j|; however, we leave it in for comparison with Nandy's calculation It is common 
to specify the forcing term through its autocorrelations 

(f Q (k, Lo)f (k', w')> = 2W(k)P ap (k)(2ir) d+1 6(k + k')5(uj + lj') , (6) 

where W{k) = Wok~ y is the forcing spectrum and the presence of the projection operator guarantees that the forcing 
is solenoidal (and hence maintains the incompressibility of the velocity field). Since the RHS is real and symmetric 
under k > — k, the configuration-space correlation is also real. 

Following FNS [2|, we impose a hard UV cut-off A <C kd, where kj is the dissipation wavenumber. With this choice of 
cut-off the theory only accounts for the largest scale behaviour (and therefore should not reproduce results for inertial 
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range turbulence). This cut-off was later relaxed to A ~ 0(kd) by YO, although the rest of their renormalization 
procedure followed FNS. The velocity held can then be decomposed into its high and low frequency modes, introducing 
a more compact notation as (u~ and it+ are often also expressed as u< and u>), 



u a {k) 



u-(k) 0<|fc|<e- £ A 
u+(k) e~ e A <\k\<A 



> 



with k = (fc,w) such that dk = (2iTY id+ ^d d k duj and 5(k) = (2ir) d+1 8(k)8(uj), allowing the NSE to be rewritten for 
the component helds: 



(7) 



(tw + v k 2 ) u a (k) = f a (k) + A M Q/37 (fe) / dj I dp \u {j)u~ (p) + 2u (j)u+(p) + u+(j)u+(p) 

8(j+p- k) 



and similarly for the high frequency modes, u~^{k). The filtered vertex operator M a ^(k) is understood to restrict 
< k < e~ e A in the non-linear term. This will later lead to an additional constraint on the loop integral. This 
constraint is neglected by many authors. 

Together with a perturbation expansion and the zero-order propagator, 



<(k) = u+ (0 Hk) + AoK+«(fc) + X 2 u+^(k) + 
1 



Go(k) = 



(8) 



iuj + i^ok 2 ' 



it is possible to solve for Ua^\k) in terms of Ua^ik) = Go(k)f^(k) using powers of Ao, which may be substituted 
back into equation ([7]). Performing a filtered-averaging procedure (•••)/, under which: 

1. Low frequency components are statistically independent of high frequency components. 

2. Low frequency components invariant under averaging: (/")/ — / _ and so (u~)f = u~ + Q(\p). 
(This can be shown more rigorously using a conditional average (22l. |23|. and discussed in [24j|.) 

3. Stirring forces are Gaussian with zero mean: (/ + )/ = 0, {/ + / + / + }/ = 0; 
and using equation ©, we obtain 



(iuj + v k 2 )u-{k)= f-(k) + Af-(k) +\ Q M-^(k) / dj dpu^(j)u-(p)6(p-k+j)+E-(k) , (9) 



where 



E-(k)=8X 2 M aP7 (k)e-(k) / dj / dp M 7 ^(p)G (p) GoQ) W(j)P^(j)u^(p + j) 



(10) 



and we have used ^-functions, 



Q-(k) = 6{Ae 



-(k) = 6{\k\-Ae- t )e{K-\k\) 



o~(p + j)o + (j)e + (p)S(p-k + j) 



1 \k\ < Ae~ e 
1/2 |fc| = Ae- e 

otherwise 

1 Ae" £ < |fc| < A 
1/2 |fc| = Ae- e 

otherwise 



(11) 



to explicitly control the shell of integration, so the momentum integrals are now < \k\ < oo. The induced random 
force, 



(12) 



A/"(fc) = XoM-^(k) / dj G (j)Mj)Mk ~ j)G (k - j)6 + (j)9 + (k - j) , 
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U a (k,(j) 

G o (k,u)P a0 (k) 



*.Q.t... = (Uk,u)f p {k>,<J)) 

k,UJ k'.Ld' 



p 



->--< 



A M Q( g 7 (fe) 



FIG. 2: Feynman rules for the velocity, force autocorrelation, propagator and vertex. These modifications clarify those in figure 
3 of 0|, since their vertex operator appears to be connected to the propagator. The dotted lines represent any combination of 
the solid lines, shown left. 



*— = — — o + 4 




FIG. 3: Renormalization of the velocity field giving the one-loop correction to the viscosity. See figure [2] for Feynman 
rules. From the factor 4 (see 2 comes from exchanging which leg from the left vertex connects to the noise 

correlation, and the other 2 from the thick line instead being incident on the LHS. 



compensates for the effect of forcing on the eliminated modes - see section|4]for more details. Note that in equation (j^J) 
we have, following Q and [j| , neglected the velocity triple non-linearity (and thus all higher non-linearities which are 
generated by it). Eyink [al showed that this operator is not irrelevant but marginal by power counting (see appendix 
however, as noted in [25| . this choice merely indicates the order of approximation and doesn't require justification. 
In any case, these higher-order operators become irrelevant as k — > Q- 

Multiplying both sides of equation ([9]) by Go (fc) and neglecting the triple non-linearity, this expression can be found 
from the graph given in figure [3] using the rules in figure [5] and the form in equation (|10[) . It should be noted that the 
symmetry factor of the graph is 4 (figure 9 in (2(| ) . 

From equation (I10[) . we may perform the frequency integrals in either order to give 

is J (2ir) d f J ,M iuj + v j 2 + v Q p z 

which, along with the definition of W(j) and M a ^ 7 (fc) = (2i) _1 P Q ^ 7 (fc), may be compared to (2.10) in [5] or (4) in 



3.1. Analysis of the self-energy integral 



We begin this section with the motivation for calling this a self-energy integral. The term has been borrowed 
from high-energy physics, and it represents the field itself modifying the potential it experiences. In high-energy 
physics, the renormalized or dressed propagator may be written using the Dyson equation (2(| (equation (27)) as 
G = Go + GqYjG + • • • , where E represents the self energy operator. In our case, we are instead writing u a (k) = 
Go(k)f a (k) + Go(k)'E af j,Uf 1 (k), where the structure S a = E^^u^ can be seen from the graph in figure [3] or equation 
([T3| once the integral over p has been trivially done. 

As can be seen in equation (|T5)) . the constraints on the integral are provided by the product of 8 + functions. We first 
show how this can be expanded before verifying that at 0(k) the substitution causes two compensating corrections, 
and hence to 0(k 2 ) there is no correction. Following this, corrections to the calculations by Wang & Wu and Nandy 
are evaluated and their contribution to the final result accounted for. 

We perform the integral over d d p so our product of 9 + functions becomes Q + (j)Q + (k — j). The second constraint, 
9 + (k — j), is sometimes ignored (see for example equation (4) in [9]) and this is a source of error in these calculations. 
With the definition + (n) = 0(\k\ — D)8(A — \n\) from equation (fTTj) where D = Ae~ e , we Taylor expand the latter 
about k — and our product becomes 

+ (j)0 + (k-j)=e+(j) (i-t3\e(A~\j\)S(\j\-D)-9(\j\-D)S(A-\j\)])+0(k 2 ) , (14) 
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FIG. 4: (Color online) The shifted shell of integration caused by the substitution used by FNS/YO. Jm and j± label the 
momenta in the directions parallel and perpendicular to that of the shift, k. (Left) The constraint 6 + {j + ife) only, as used by 

Wang & Wu; (right) + (j + \k)6 + (j - ±fc) as used by FNS. White ( ) shows the position of the unshifted shell 6 + (j), with 

light grey shading ( — ) showing the shifted shell(s). Small tick marks indicate the centre for the shifted shell(s) at j = ^fhk. 
The dark grey highlights the overlap of two shifted shells. The green point (in the top-left quadrant) is a random momentum, 
q, which lies within the resultant shell, while the red cross (lower-right quadrant) shows — q. Clearly, in the left case the shell 
is not symmetric under j — > —j, and as such identities requiring a symmetric shell are invalid. By correctly accounting for 
the additional momentum constraint, we are led instead to a shell like the one to the right, which is symmetric. In the right 
figure, where ||fe = is shown by the (horizontal) height of the orange-filled triangles towards the left, we see that as k — > 
(the triangles shrink to their vertical baselines) the overlap increases and the resultant shifted shell is well approximated by 
the original, unshifted shell (• — •), as found in equation (|17|l . 



see appcndix|B]for details. We see that the additional constraint has introduced a first order correction to the constraint 
on \j\. Further, the presence of the (5-functions show that these contributions are evaluated on the boundaries. This 
correction is absent from the work of Wang & Wu, Teodorovich and Nandy as they ignore this constraint. We shall 
see later that, from a diagrammatic point of view, this is equivalent to ensuring that all internal lines have momenta 
in the eliminated band. 



3.1.1. FNS and YO 



We now turn our attention to the substitution j — > + j made by FNS and YO, under which our constraints 
become 



9+(j)9+(k-j) ^0+(\k + j)9+{\k-j) . 
Taylor expansion of these high-pass filters is now 



ak±j) = e+(j)± 



kj 



2j 



0(A - j)S(j — D) — 9(j - D)S(A - j) + Oik 2 ) 



-(j)±x(k,j) + 0(k 2 ) , 



and the product becomes 



6+{\k 



-(ifc 



-(j) + 0(k 2 ) 



(15) 



(16) 



(17) 



The contributions at 0(k) cancel one another exactly, and there is no correction to the simple constraint on 
Without the constraint 9 + (k — j), the substitution would have led to just 9 + (^k + j), which clearly does introduce 
a first order correction. These points can be seen in figure |U Using this, we go on to find the result of Yakhot & 
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Orszag (a generalization of the FNS result) in the limit ui — > 0, k — > 

» / \ S d d 2 -d-eW Q X 2 fe ee -l 

A^ (0,0) - 



(2ir) d 2d(d + 2) v 2 fc 
^ d A 2 (0) f— , (18) 



where 

1 ~ (2^^ ' ™ a ~ 2d{d + 2) ' 



Ad = 7K ^ S A d , A d = ———, A 2 (0) = ^-^, e=4-d+y. (19) 



3.1.2. Wang & Wu 



Wang & Wu were unsatisfied with the substitution used by FNS and YO. This is because they do not impose 
the condition that D < \k — j\ < A on the self-energy integral, and on the face of things the substitution shifts the 
integration domain (see figured]). The authors then continue without making any substitution but simply Taylor 
expanding the denominator and expanding the vertex operator to 0{k 2 ) 



4A 2 

-^-M aPl (k)u-(k,0) 



= *Wfc) / d d j r y - 4 e + (j)9 + (k-j) 



k v P^(j)P 0v (j) - kJ-^P^(j) - k v J -^P M {j) - j„P fi jU) 



3 i±3v 



(20) 



where the operator P a g~„{k) 



W Xl M aM (k)u^ (fe,Q) 



is defined for convenience. We now expand the product of 9 + 



functions as in equation (fT4")) but note that only the last term in the square brackets above is not already of 0(k). As 
such, it is the only term that can generate a correction. £~ then splits into 



6Z- 



where 

is the contribution used by Wang & Wu without imposing the additional constraint, and the correction 



WO ld d jj- y ~ 4 e+(j) 



= / dd J r y - A 0+(j) J -^P Pl (j) [0(A - j)S(j — D) — 9(j D)S(A j) 

includes the additional boundary terms. 

The first contribution above leads to the Wang & Wu result, that 



where 



A^«-(0,0) = i/ ^(0) 



-Av^ w (k,0)k 2 u-{k,0) , 



A* - Sd A* 
Ad ~ (2n) d d 



d 2(d + 2) 



(21) 



(22) 



(23) 



(24) 



(25) 



We now evaluate the first order correction given by equation 
(see, for example, (27l.[28j). 



p 



Sd 



; v d(d + 2) 

WoAj 1_ 

vl (27r) d 2d(d + 2) V A 



(d + 1)6^6^ — S^pfivy — S^Svp 



using the standard convention that 9(0) = | 
D~ e - A _e l 0(0) 



(26) 
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and so the correction to the viscosity increment found by Wang & Wu is 

5v ww (k o) - W ^ g? 1 (flzT\ (27) 
^° (fe '° j - ^ 2 (2^2d(d + 2) ^ A* J' (27j 

If this contribution is added to the result for the renormalized viscosity increment found by Wang & Wu in equation 
([25]). we find 



Avq(0, 0) = Av™ w (0, 0) + 5vff w (0, 0) 
S A ( d-1 1 



° A °fe rf -l) 



(27r) d V2(rf + 2)e 2d(d + 2) / v ; (28) 

_ S d d 2 ~d-eW \ 2 ( e^- V 
~ (27r) d 2d(d + 2) z/ 2 A £ ^ T 

which is exactly the result obtained by YO, see equation (fT8|) . Hence we have shown that a more careful consideration 
of the region of integration used by Wang & Wu instead leads to the result found by FNS and then later YO. The 
approach taken by Teodorovich [6[ uses a different method for evaluating the angular part of the self-energy integral. 
However, the author misses the same constraint and thus arrives at the result as Wang & Wu. 



3.1.3. Nandy 

In the paper by Nandy [111 ], the author presents an argument based on symmetrising the self-energy integral. 
Referring to equation (TT5|) . he points out that there is no reason to do the d d p integral first, and that the result should 
be an average of the two. Performing the d d j integrals first, gives 



A\2 

E-(fc,0) = -£Map y (k)u-(k,0) 

v 



^(l*-|>M_(p)#^^- P >e+M 



(27r) d \k-p\- 



\k — p\ 2 + p 2 



(29) 



Taylor expanding the function \k — p\ y 2 and the denominator, and using the definition and properties of the vertex 
and projection operators, to 0(k 2 ) leads to 



= *W fe ) / d d PP - y - 4 9 + (k-p)9 + (p) 



k u P lv {pf-^ - k v P^{p) V -^L + kpPw(p) + (y + i)KP M {p) P -^+p,P Pl {p) 



(30) 



Again, we see that all the terms in the square brackets are already O(k) except the last one, and so this is the only 
term which generates a correction. Once again decomposing 



(31) 



the contribution calculated by Nandy is 



kp^pP^- - kp^M^t- + W/» + (y + i)k v p til {p) p -^+pM P ) 



(32) 



and the correction generated by expanding the product of 9 + functions is given by 



= -PafruWK / d d pp 



{ P Y^Lp M {p) 9(A - p)8(p -D)- 9{p - D)6(A - p) 



= -Si:; 



(33) 



We see that, with the relabelling p —> j, the correction is exactly the same as equation (1231) only with the opposite 
sign. Therefore, we see that 



i[< Wang& Wu> + <Nandy >] = i [£~ + £~] = \ [(£" + ST.-) + (£" + ST')] = T' c 



(34) 
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which is why this symmetrization produced the correct result. 
Nandy for performing the integrals in this order 



In fact, evaluation of E a leads to the result found by 



A^(0,0) 



WnX 2 S d 



d-2e 



v 2 {2jr) d 2d{d + 2) 



1 



d 2 - d - 2e = d 2 + d - 2y - 8 



and the correction is 



H"(o,o) 



W>A 2 S d 



{2n) d 2d(d + 2) V eA 



(35) 



(36) 



Combining these results we again find the result of Yakhot & Orszag, equation (fT8"|). showing that regardless of which 
integral is performed first we obtain the same result and so a symmetrization is not necessary. 

In a completely different approach, Sukoriansky et al. 24) used a self-substitution method to solve for the low- 



frequency modes and claim to evaluate the cross-term exactly. This method does not generate the cubic linearity, 
instead it creates a contribution of the same form as Wang & Wu equation (I2TH) but with the condition 0~(k — j). 
Combined, this then covers the whole domain, and the authors drop any conditions on |fe — j\. However, 



0- (k - j) + 9+(k- j) = 0{k - \k - j\) + 1 , 



(37) 



due to the upper momentum cutoff. This can then be Taylor expanded for small k, and leads to a contribution from 
the upper boundary, neglected in their analysis. In fact, this correction places their result somewhat between that 
found by Yakhot & Orszag and Wang & Wu, 



A^o(O,0) 



Wo A 2 



Wo A 2 



A* 



(Ae" 



A- 



2d{d + 2) 



A*(Ae- e yc - A d (e)A~ 



(38) 



where it is missing the contribution from the lower boundary. However, this self-substitution is not the same as 
solving a dynamical equation for the low-frequency modes and substituting for the high-frequency components. This 
method is fundamentally different from the standard RG procedure, and its result agreeing with neither YO or WW 
is further evidence that it is another approximation entirely. 



4. NOISE RENORMALIZATION 



4.1. FNS treatment 

In the paper by FNS [2J , the authors use two different scaling conditions (see appendix [Aj when analysing their 
models due to the contribution of the induced force to the renormalization. We first mention the results used by 
Yakhot & Orszag for comparison (see appendix |A"|) . 

X= ±(3z + d + y) , x = d+l — > * = 2-§, X* 2 = , (39) 

where A* 2 is the reduced coupling at the non-trivial fixed point. 

FNS model A (y = —2): For this model the authors show using diagrams (see figure 1 of Q) how the propagator, 
force autocorrelation (shown here in figure [6J and vertex are renormalized. They conclude that vq and Wo are 
renormalized the same way in their equations (3.10-11). This condition implies fixing the mean dissipation rate 
rather than Wo , and is then enforced under rescaling (see appendix |Aj by choosing \ = z + which does not agree 
with x above (first relation) used by YO with y = — 2. 

At this point, FNS invoke Galilean invariance (GI) to impose the condition that the vertex is not renormalized, 
such that A — Ao = 1 to all orders in perturbation theory (in the limit of small external momenta, appendix B of 
) . While this is the case at k = j29| , and as such does not invalidate the FNS theory, in general the consequences 
of the symmetry are trivial and do not lead to a condition on the vertex [291 - 1321 ] . Further discussion is given in the 
conclusions, Sec. [5j Taking the condition x = d + 1 to preserve Galilean invariance, they find the non-trivial stable 
fixed point A* 2 = e/2A d (when e > 0). 

FNS model B (y — 0): In this case, the one-loop graph in figureHlis claimed to be 0(k 2 ) and so cannot contribute 
to the constant part of the force autocorrelation. This term is then irrelevant and the force is rescaled accordingly. 
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AoAVv(fe) / dqG (q)f+(q)f+(k - q)G (k - q) 



FIG. 5: Graphical representation of the form of the induced random force given in equation (|12|l . Note that the two momenta 
q, k — q lie in the eliminated shell. 




= < </* f$) 4 + 2u k a A aP (k,k') U $ 
FIG. 6: Renormalization of the force autocorrelation by FNS. See figure[2]for Feynman rules. 



This requires \ = \{'& z + d), which is the same condition found by YO with y = 0. Ensuring that Galilean invariance 
is satisfied, they have A* 2 = e/3A<i, as do YO. 

This difference in scaling conditions leads to different differential equations and hence different solutions for the 
reduced coupling A and the viscosity, depending on whether the noise is allowed to be renormalized or not. In the 
field-theoretic approach by Ronis [lj|, the force is also allowed to be renormalized, and the author comments that 
YO ignore this in their analysis. In fact, they restricted their work to y > —2 to avoid this issue. 

This discrepancy only really applies to y = —2 when the noise coefficient is renormalized, although could lead to 
complications for y < —2 as the induced force always contributes as 0(k 2 ) to the autocorrelation and becomes the 
leading order as k — > 0. In their paper YO state that "in the limit k — > this [induced] force is negligible in 
comparison with original forcing with y < —2", and present an argument for neglecting it as equation (3.13). For the 
case y > —2 it is sub-leading and thus safely neglected. This highlights another potential problem with calculating 
inertial range statistics, since it is only sub-leading as k — > 0. 



4.2. Re-evaluation 



In this section we will show how the induced force leads exactly to the graph in figure HJ We then evaluate the 
graph to analyse the contribution to the renormalization of the force. For this, consider the form of the induced force 
shown in figure [5] Under averaging, we see that the graph forms a closed loop and is 0{k) due to the vertex operator, 
and hence (A/ Q (fc)) = 0. The new random force f„(k) = fa(k) + is invariant under the filtered-averaging 

procedure, and has autocorrelation 

(f a (k)f P (k')) = (fa(k)f P (k')) + (Af a (k)Af p (k')) • (40) 
The new contribution due to the induced force is written as equation (|4T[) . 



(Af a (k)Af p (k')) = \lM aM {k)M p ^{k') J J dqdp G (q)G (k - q)G (p)G (k' -p) 

x (Mq)Mk-q)U(p)f v (k' -p))6+(q)6+(p)6+(k-q)9+(k' - p) . 
Since the forcing is taken to be Gaussian, we may split the fourth-order moment 

and, using the definition of the force correlation equation ([SJ, we see that the first contribution leads to two discon- 
nected loops (which do not contribute to the force renormalization), whereas the other two both generate graphs like 
figure [6] and appear to contribute towards the renormalization of Wq. 
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Using the rules given in figure [2 we write an analytic form for the C(A§) diagram in figure [6] as equation (|43|) . 



A ap (k,k') = 2\ 2 M a01 (k)M pp „(k') J J dqdp Go(q)(fp{q).Up))Go(p) 

Go(k - q)(Mk - q)f v (k' -p))G (k' -p) . 



(43) 



The factor 2 is due to symmetry of the graph (figure 2 in [26j]). This may be compared to the correlation of the 
induced force given by equation ()41[) which, along with the requirement that momenta of all internal lines in equation 
(|4"3"|) are in the eliminated shell, agree exactly. 

An outline of the evaluation of this correction to leading order is given in appendix [Cj As a result of our analysis, 
wc find 



(fa(k)f P (k')) = 2W k-yp ap (k)5(k + k') 



i + A 2 (o)s d [ —L i ky +2 



(e + y + 2)Al 



y+2 



(44) 



where 



B d 



S d d 2 -2 
(2n) d 2d{d + 2) 



(45) 



We see for the equilibrium case y = —2 that the correction may be taken as an multiplicative renormalization to 
Wo and write 



(fa(k)Uk')) = 2W z k 2 P ap (k)S(k + k') , 



(46) 



with 



Wi = W 



l + \ 2 (0)B d 



(47) 



This may be compared to equation (3.11) of FNS [2j and (3.4b) of Ronis [l8[ (with y — —2). We note that both 
authors find the noise coefficient and viscosity to be renormalized with the same prefactor, what we have defined as 
B d - However, this prefactor only coincides with A d found by FNS and YO for y = —2 (e = 2 — d), and the analysis 
is only valid for the equilibrium case (the work by Ronis is an expansion about y = —2). 

The prefactor B d was calculated to leading order with no change of variables in a similar fashion to Wang & Wu, 
and agrees with A d found by FNS and that by YO with y = — 2. If the prefactor A* d found by Wang & Wu and others, 
which is e-independent, were the true expression, we should have recovered it from this analysis also. Instead, it only 
agrees with B d when d = 2 (the critical dimension for y = — 2 (where e = and A * d = A d )). Taking the induced 
contribution as an multiplicative renormalization only when y = — 2, A d = B d whereas only A* d=2 — B d= 2- We feel 
this supports our argument that the correct expression for the prefactor is the e-dependent result found by YO. 

As noted by Ronis, the pole structure leading to logarithmic divergence in the noise renormalization more generally 
occurs when e + y + 2 = 0. However for the case y = —2, we recover the same pole in e found by FNS and in the 
viscosity renormalization. Since noise renormalization is only meaningful in this case, the presentation of equation 
(3.4b) in [l~8[ as a general result seems misleading. 

In summary, we have calculated the renormalized noise coefficient to one-loop and find the prefactor B d to agree 
with FNS as well as YO when setting y = —2 in their result. While noise renormalization was not considered by 
Wang & Wu, we have found their e-free result, A d , to only agree with B d in 2-dimensions. For y > —2, this induced 
forced correlation becomes sub-leading and is ignored, as assumed in the YO analysis. When y < —2, the induced 
contribution does not renormalize the noise coefficient but will be the leading term as k — > 0. In this case, it is not 
clear how to interpret the validity of the results obtained, since the forcing appears on large scales to be dominated 
by the order A 2 (0) contribution, making the viscosity calculation order A 4 (0), i.e. two-loop, which has not been done 
here. 



5. CONCLUSION 



A summary of our results and a comparison with other authors is presented in table HI We conclude that the 
analysis of FNS does not suffer from a shifted domain of integration in the self-energy integral which is evaluated 
due to the constraint \k — j\ < A, neglected by other authors. Using 9 functions to control the integration domain, 
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Our analysis 


FNS A 


FNS B 


YO 


WW / T / N 


R,onis 


Viscosity 


id(e) 




a \ >ly = 


A d (e) 


id = Ad(0) 


Ad(e)\ 

"V /\y—_2 


Noise 








B d 




B d 


Pole structure 


l/(e + y + 2) 


1/e 




iA 




l/(e + y + 2) 



TABLE I: A summary of the prefactors for the viscosity increment, noise renormalization, and the e-pole structure found in the 
various analyses considered in this paper. These expressions are valid for all y, with the exception of FNS models A (y = —2) 
and B (y = 0). T represents Teodorovich and 'N' Nandy. A d (e) = (d 2 - d - e)/2d(d + 2) and hence B d = (d 2 - 2)/2d(d + 2). 
By pole structure we mean the e-dependence of the denominator for the induced noise correlation (see equations (|18[) and (144 1) , 
for example). Our analysis agrees with FNS and YO. The work of Ronis appears to be for general y, but the viscosity is only 
in agreement for y — — 2. 

we have shown that the corrections cancel exactly at first order in k when the change of variables is made. We then 
showed that this ignored constraint leads to a correction in the Wang & Wu- and Nandy-style calculations which 
exactly reproduces the result found by YO. The noise renormalization for the case y = — 2 was then shown, using 
a substitution-free method similar to Wang & Wu, to lead to a prefactor compatible with YO for all d and only 
compatible with Wang & Wu for d = 2, which we feel supports our claim as to the validity of the FNS and YO results. 

That said, some comments should be made on the application of this method to calculating inertial range statistics, 
which may not be so well justified. Despite its applicability only on the largest of length scales, Yakhot & Orszag use 
the expressions obtained with this infra-red procedure to calculate inertial range properties 5] , such as the Kolmogorov 
constant. To do this, they use a set of assumptions that they term the correspondence principle. 

Briefly, the correspondence principle states that an unforced system which started from some initial conditions with 
a developed inertial range is statistically equivalent to a system forced in such a way as to generate the same scaling 
exponents. In particular if forcing is introduced to generate the scaling exponents at low k, this artificially generated 
"inertial range" can then be used to calculate values for various inertial range parameters using the properties of 
universality. There is an implicit assumption that, as long as the scaling exponents match, all other quantities will 
also match. This may be the reason that YO raise the cutoff A out of the production range to O(kd) (see above 
equation (2.2) in so that the renormalization passes through the inertial range, whereas FNS explicitly consider 
A <C kd (final paragraph of their section II. A). 

YO find that when y = d the noise coefficient, Wo, has the dimensions of the dissipation rate, e, and they take 
Wq = ae (with a constant). They can then obtain a Kolmogorov scaling region when e = 4 is used, but also require 
e = in the prefactor Ad{e) in the same equation. This has been unsatisfactory for many authors, and appears to 
favour the e-free result found by Wang & Wu, as then e = 4 alone reproduces the famous fc -5 / 3 result. However, we 
have shown why the e-free result is incorrect. 

There are still a number of technical difficulties associated with taking e — > 4 and generating a fc -5 / 3 spectrum: 

• The Wilson-style e expansion is valid only for e small, and there is no evidence that results will be valid at 
e = 4. The neglected cubic and higher-order non-linear terms generated by iterating this procedure may not be 
irrelevant, and there is no estimate of the accumulation of error even for e ~ 0, let alone e — > 4. In the review by 
Smith & Woodruff [33[ , they discuss the only justification for the validity of e — > 4 being that it leads to good 
agreement with inertial range constants, and describe it as "intriguing and difficult to interpret" . They also 
present an argument for YO's use of e = in the prefactor, it being required for a self-consistent asymptotic 
expansion at each iteration step. 

• The IR behaviour as k — > is dominated by the fixed point which, for e > 0, is at A* = {e/iAd) 1 ^ 2 ■ To lowest 
order in e, this is then evaluated with Ad(e = 0). However, e is no longer small, nor is A*. In 3-dimensions with 
e = 4, this fixed point is at A* = (20"7r 2 e/(6 — e)) 1 ^ 2 — 11.5 to leading order in e, or ~ 20 when evaluated to all 
orders. 

• As shown by figure [T] the asymptotic nature of this renormalization scheme taking us to the infra-red means we 
don't enter the inertial range, and are always sensitive to the forcing spectrum. 

• The forcing spectrum required to obtain fc~ 5 / 3 is divergent as k — > (e = 4 requires y = d > 1, so W(k) ~ k~ d ), 
as is the energy spectrum itself. As shown by McComb |7[ , ensuring that there is a balance between energy input 
at large length scales and energy dissipated at small (this is statistically stationary turbulence) we see that the 
range of forced wavenumbers predicted by their analysis has kt/kb — 1.007, where k t and kb are, respectively, the 
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upper and lower bounding wavenumbers of the input range. The energy input is also logarithmically divergent 
as k t — > oo or fcj — > 0. 

• The condition of Galilean invariance (GI) used by FNS and adopted by YO to enforce the non-renormalization 
of the vertex at all orders is actually only valid at k = [29| . In general, the consequences of GI arc trivial 
and provide no constraint on the vertex [29l - l32j . This is supported by recent numerical results [34| from a KPZ 
model on a discretized lattice with a broken GI symmetry, which have found the same critical exponents as the 
actual KPZ model (which does possess GI) @, Il4j . even though GI has been explicitly violated. This questions 
the connection between GI and the scaling relations associated with the critical exponents. As such, care must 
be taken when extending this theory to k ^ 0. This introduces another issue for the study of inertial range 
properties using the correspondence principle, as k cannot be chosen to lie in the inertial range without the 
vertex being renormalized. 

• The assumed Gaussian lowest-order behaviour of the fluid is only valid at the smallest wavenumbers when 
subject to Gaussian forcing, since the response of the system is then also Gaussian. However, this assumption 
cannot be translated to the inertial range, which should be insensitive to the details of the energy input and is 
inherently non-Gaussian Q. 

The need to use two different values for e in the same formula to estimate inertial range properties is therefore not 
the only failure of this scheme. 

The solution for the renormalized viscosity at the largest scales (I — > oo) can be found to behave as 



where // = Ae — > is the new cut off. With the assumption Wo = as, interestingly this does have the same form as 
that found by other methods (e.g. Q), that the viscosity is proportional to e 1 / 3 and, with e = 4, the cut off /i -4 / 3 . 
However, there is an important difference: the cut off /i is going to zero for this expression to hold, which is not the 
location of the inertial range, unlike the iterative averaging approach by McComb QjlUlj]. Smith & Woodruff note 
[33l ] that this is not dependent on the dissipation range quantities A and vq, as inertial range coefficients should be. 
But it is still dependent on the forcing spectrum through e, which it should not be. 



We thank David McComb for suggesting this project and for his helpful advice during it. We were funded by STFC. 



Although rescaling the variables after performing an iteration of the renormalization procedure outlined above is 
not performed in the calculation of the renormalized viscosity (and thus it can be argued not to be an RG procedure 
Q), it is still useful to consider how the rescaling would affect the equations of motion. Using a scaling factor, s, 
the spatial coordinates transform as x — sx' and t — s z t' (where the unprimed variables are the original scale), with 



s > 1, and so k = s^k', lo = s- z uj' with u(k,u) = s*u'(k',u'). In YO [|, s = e £ , a(£) = zl and ({£) = = S*. 
Equation ([3|) then transforms under the scaling to 




(48) 
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Appendix A: Operator scaling 




(Al) 



and so we find 



v' = v{s) 
V = \{s) 




(A2) 
(A3) 



_ s z -x+i<-y+ d + z )f a (k',u J ') . 



(A4) 
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Using equation (|A4[) with the definition of the force autocorrelations equation © and the scaling for k and u>, we find 

W' = W(s) = s 3z-lX+d+y w ( A5 ) 



Equations (|A2[) — (|A5|) agree with (2.28)-(2.33) of |5|. Due to Galilean invariance, as k — > equation (|A3[) is forced 
to give the condition that x = d + 1. For y ^ — 2, the elimination of scales should not affect Wo as there is no 
multiplicative renormalization, a condition that must also be preserved under scaling to find \ — \ (3z + d + y). YO 
note (from (2.34) of [f|) that the renormalized viscosity at the fixed point is s-independent if z = 2 — AjX* 2 = 2 — |. 

As noted by FNS and discussed in section [U for the case y = —2 we must consider the renormalization of Wq and 
instead require that v and W be renormalized in the same way, i.e. x — z + i- This scaling condition is not the same 
as that for y ^ —2, and leads to a different solution. 

Under the YO prescription, the triple non-linearity 

(A6) 



liM aM {k) I dj / dp M lllv {k-j)G Q (k-j)up{j)u^{p)u u {k-j-p) 



gives, using the expression \ — \ (3z + d + y), 



This should be compared to equation (2.45) in Q, which reads 



(A7) 



(A8) 



They comment that for y < d (e < 4) the operator is irrelevant, and marginal when y = d (e = 4). However, we see 
that their result requires x = h{d + 2/ + 2), which only agrees with the above expression for x (ensuring Wq is not 
altered) when z = | = (2 — |)| . Therefore, they have already used e = 4 to obtain this result. If we do not specify 



e but do require that z = 2 — | (so that the viscosity at the fixed point is s-independent), we see from equation (|A7[) 
that // ~ /i and the operator is not irrelevant but marginal. (This could also have been seen by requiring x — d + 1 
in equation (|A7[) . and we see that if the vertex is not renormalized the triple moment cannot be irrelevant.) This is 
discussed in a pape r by Eyink Q . Attempts to retain the effects of the triple non-linearity on the viscosity increment 
are analysed in [35l - l37j . 



Appendix B: Taylor expansion of # + -functions 



We here describe the procedure for Taylor expanding a # + -function. The high-band filter 9 + is defined as 

+ (j)=6(\j\-A)9(A-\j\) , A = Ae- l <A, (Bl) 
where the first restricts us to j > A and the second to j < A. The ^-function product of consideration here is 

9+(j)9+(k-j) , (B2) 
and we Taylor expand 9 + (k — j) = 9(\k — j\ — A)9(A — \k — j\) as: 

:• (V9(\k-j\ - A)) +... 

k=0 



\k-j\-A) = 9(\k-j\-A) +k-[V9{\k-j\-A) 

k=0 \ 



= 0Qj\-A) + k- 



\k-j\ 



6(\k-j\-A) 



k=0 



k- j 



9(\j\-A)-—fS(\j\-A) + 0(k') 



(B3) 



9(A -\k- j\) = 9(A -\k- j\) + k ■ (V9(A -\k- 

k—0 \ / 

6{A-\j\)-k- ( k ~ 3 



\k-j\ 



(s(A-\k-j\) 



k=0 



k=0 



0(A-\j\) 



kj 



S(A-\j\) + 0(k 2 ) 



(B4) 
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Our expansion is then 



e + (k-j) = e + (j) 



kj 



9(A - j)S(J -A)- 6{j - A)S(A - j) + 0(k 2 ) . 



(B5) 



Appendix C: Evaluation of the noise renormalization 

We now evaluate the correlation of the induced force to leading order using a more compact notation. Starting 
from equations (PTT1H!?)) . 



(CI) 



'q u p u k—q k f —p 



we note that our integrals here are unconstrained and the shell of integration is controlled by the # + -functions. Since 
the substitution p — »■ p' = k' — p preserves the product 9 + {p)9 + (k' — p) — > 9 + (p')9 + (k 1 — p') = 9 + {k' — p)9 + (p), we 
may use it, along with the property of the vertex operator M pliV (k') — M pv ^(k') and index relabelling for the second 
term in the square brackets, to combine the two contributions and write 



<A/*A/*') = 2X 2 M^M^ / dqdp GlG^GtG^v{ffji){f k 1 -«ft'- p )9+9 p - 



+a+ a+ 

k—q k' —p 



(C2) 



which we see is exactly equation (|43|) . This reveals that the symmetry factor of 2 associated to the graph is due to 
exchanging legs on the vertex. In to this we substitute the definition of the force autocorrelation, then integration 
over p is trivially done using the (5-functions obtained to give 



<A/ Q fc A/ p fc '} = 8\lM^ 7 M«J(k + k') / da 



G, 



k—q 



w( q )w(\k - q\)P^p* v - q e+e+_ q . 



(C3) 



The constraint enforced by the remaining 8- function is then used to restrict k' = —k, along with the property 
M pfiU {k')8(k + k') = ~M ppv {k)5{k + k'), resulting in 



<A/ Q fc A/ p fc '> = -8X 2 Q M^M^J(k + k') j da 



Gl 



G k -t 



w( q )w(\k - q\)P q ^P*- q e+e+_ q . 



(CM) 



The frequency integral is then performed, closing the contour in the upper-halfplane and collecting the residue from 
two poles ft = iv^q 2 and Q = u + ivo\k — q\, with the result 



/ 



dCl 



Lr 



G 



k—q 



q 2 (iu> + u Q (q 2 + \k- q\ 2 )) + |fc - q\ 2 (iu - ^{g 2 + |fc - q\ 2 )) 
q2\k-q\ 2 (ij 2 + v$(q 2 + \k-q\ 2 ) 2 ) 
1 



iuj + V{){q 2 — \k - q\ 2 
1 



v '6 iq 2 \k - q\ 2 (l 2 + \k — q\ 2 ) 

The limit to — > offers a huge simplification to the result. This is inserted in to equation (|C4 

id 



<A/*A/f) = ^^M^M^J(k + k')J (q\k q \yy- 2 



|fc_„|2 q k-q 



and the integrand is expanded to leading order in k as k — >• 

\2 W 2 r d d q 



(A/ Q fe Af fe ') = 



} = JW^^ + y } J -^L d q- 2 ^)P^ V 9+ + £>(*») ■ 



(C5) 



(C6) 



(C7) 



Note that there is a power of k associated to each of the vertex operators, hence the leading contribution will always 
go as k 2 . Expanding the function 9 + {k — q) we do not generate corrections as we are working to zero-order in k in 
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the integrand and the corrections are 0(k). Expanding the projection operators and performing the (d — 1) angular 
integrals we find 



J dq q- 2 ^+ d - 1 9+ , 

where we expand the vertex operators, do the remaining integral and perform contractions to obtain 



which we rearrange to our final result 



J{e+y+2) _ 1 \ ^ _ f f( _ g d d 2_ 2 



(C8) 



\ 2 W 2 S* / P £ ( e +f+ 2 ) — 1 \ 

(A/ a ( fe ,0)A/^0)) = ^^ + ^^^[2^( fe )(^-2)] L £ + y + 2)K+y+ . ) ■ (C9) 



(A/ Q (fc, 0)A/ p (fc', 0)) = 2W \ 2 (0)B d [ (e + y + 2)A y +2 1 k 2 P ap (k)S(k + fc') , 5, = ^ — - - . (CIO) 
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